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ABSTRACT 


The  raagnetocaloric  refrigerator  uses  the  adiabatic 
temperature  response  of  ferromagnetic  materials  (rare- 
earth  metals)  to  a  high  magnetic  field  as  a  basis  for 
operation.  A  mathematical  model  of  a  reciprocating 
magnetic  refrigerator  having  an  active  element  of  a 
rare-earth  metal  predicted  that  a  60-K  temperature  span 
could  be  achieved.  Based  on  the  model's  results,  an 
experimental  refrigerator  was  fabricated  and  the  thermal 
characteristics  of  an  active  regenerator  in  a  recipro¬ 
cating  apparatus  were  measured.  The  active  material 
was  gadolinium,  which  has  a  Curie  temperature  of  293  K 
and  was  formed  into  an  embossed  ribbon  geometry.  Signi¬ 
ficant  regeneration  in  the  magnetic  material  was 
observed,  and  a  50-K  temperature  span  was  received.  A 
cold  end  temperature  of  246  K  was  measured. 

The  results  indicated  that  the  active  regeneration 
concept  can  successfully  produce  magnetocaloric  refrig¬ 
eration  at  temperatures  near  the  Curie  temperature  of 
the  material,  and  the  measured  temperature  span  verified 
the  mathematical  model  developed  to  describe  the  pro¬ 
cess  . 

Further,  if  the  active  magnetic  element  were  a 
combination  of  several  rare-earth  metals  having  signifi¬ 
cant  magnetocaloric  effects  at  decreasing  temperatures , 
the  possibility  exists  for  producing  a  single-stage  mag¬ 
netic  refrigerator  that  could  operate  over  a  very  large 
temperature  span. 

ADMINISTRATIVE  INFORMATION 

This  work  was  performed  under  the  David  Taylor  Naval  Ship  Research  end  Devel¬ 
opment  Center's  IR/IED  project  entitled,  "Magnetocaloric  Cryogenic  Refrigeration," 
Work  Unit  2712-150,  Program  Element  61152N,  and  Task  Area  2R0240101. 

INTRODUCTION 

The  Navy  at  the  David  Taylor  Naval  Ship  Research  and  Development  Center 
(DTNSRDC)  has  in  the  past  twelve  years  been  developing  the  neoess.v  y  machine rv  { or 
use  with  advanced  electric  propulsion  systems  for  Suture  ship-..  Ihi-  c>ou  i  pment  use 
the  advanced  technology  ut  superconductivity  i>>  improve  tin  power  dec.  si  tv  ut  large 
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trie  propulsion  motors.  Superconductivity  occurs  when  certain  materials  are  cooled 


to  very  low  temperatures  (typically  5  to  20  K)  and  electrical  resistance  drops  to 
zero.  This  phenomenon  has  been  demonstrated  to  be  very  valuable  for  building 
highly  efficient,  high  power  density  electrical  machinery.  Superconducting  motors 
and  generators  require  a  cryogenic  refrigerator  to  maintain  the  superconductor  at 
these  low  temperatures.  Cryogenic  refrigerators  have  been  used  extensively  in  the 
laboratory  where  efficiency  and  reliability  are  of  secondary  importance.  For  naval 
shipboard  applications  the  unit  must  be  efficient,  rugged,  and  reliable  to  meet 
the  Navy's  operational  requirements.  Areas  of  concern  associated  with  these  refrig¬ 
eration  systems  are  the  compressor  inefficiency,  oil  contamination,  and  engine 
valves  and  seals  in  the  liquifier.  Thus,  having  a  shipboard  liquifier  that  required 
no  compressors  would  be  highly  desirable,  not  only  to  eliminate  the  inefficiency 
of  the  compressors  but  to  improve  the  operating  conditions  to  which  the  refrigerant 
is  exposed.  The  compact  and  rugged  construction  of  a  magnetic  refrigerator  would 
also  allow  for  operating  under  high  shock  and  vibration  loading  such  as  that 
encountered  in  a  Navy  shipboard  environment. 

Until  recently,  the  magnetic  refrigerator  had  been  used  only  at  temperatures 
below  4  K.  But  the  concept  of  magnetic  refrigeration  has  shown  some  real  possibil¬ 
ity  at  temperatures  greater  than  4  K  with  the  use  of  ferromagnetic  material.  The 
magnetic  refrigerator  uses  the  magnetocaloric  effect,  which  allows  certain  materials 
to  undergo  a  large  adiabatic  temperature  change  or  isothermal  entropy  change  when 
subjected  to  a  large  magnetic  field.  Using  this  effect,  it  has  been  proposed  that  a 
refrigerator  could  be  designed  to  have  a  Carnot  efficiency  of  about  50  to  70% 
[Barclay,  1].  In  fact,  Lacaze  and  his  co-workers  [2]  measured  efficiencies  of  64% 
in  a  double-acting  magnetic  refrigerator  at  2  K.  This  compares  to  efficiencies  of 
about  1  to  10%  for  gas  cycle  refrigerators  near  liquid  helium  temperatures  (4.2  K) 
as  reported  by  Strobridge  [3]. 
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Magnetic  refrigerators  can  be  classified  into  four  categories:  (1)  continuous 
flow  (rotary);  (2)  periodic  flow  (reciprocating);  (3)  nonregenerative  (small  AT); 
and  (4)  regenerative  (large  AT).  The  rotary  concept  (continuous  flow)  is  being 
investigated  by  Barclay  [4]  at  Los  Alamos  National  Laboratory.  Dr.  Barclay's  work 
is  in  the  4  to  20  K  temperature  range  with  no  regeneration.  The  reciprocating  con¬ 
cept  was  investigated  by  Brown  j5]  at  the  National  Aeronautics  and  Space  Administra¬ 
tion's  Lewis  Laboratory.  He  used  a  regeneration  system  in  the  fluid  to  demonstrate 
a  large  (80-K)  temperature  span.  This  worx  was  performed  near  room  temperature  to 
make  use  of  fluids  with  high  heat  capacities.  A  second  method  to  obtain  regeneration 
is  to  do  it  in  the  active  ferromagnetic,  or  paramagnetic  material,  a  method  described 
as  an  active  regenerator  in  a  patent  application  by  Barclay  [A], 

Because  the  magnetic  refrigeration  process  is  reversible,  the  thermal  and  mag¬ 
netic  characteristics  of  the  materials  vised  will  determine  how  well  the  concept 
works  and  thus  determine  the  efficiency  of  the  cycle.  The  key  to  the  process  is 
thorough  understanding  of  the  thermal  and  magnetic  interactions,  so  as  to  keep 
entropy  generation  to  a  minimum.  At  DTNSKDC,  we  decided  to  confine  the  thermal  and 
magnetic  interactions  to  one  location  by  use  of  an  active  regnerator  in  a  recipro¬ 
cating  configuration.  A  mathematical  model  of  such,  a  svstem  was  formulated,  a 
prototype  system  was  built  and  tested,  and  results  oi  the  model  were  compared  with 
expe  r i ment a  1  measu re ment s . 

BASIC  PR  I  NO  l  PLUS  OK  MAOV-iTiC  KP  i-  ,  IKK  A  lit  *. 

!'he  magnetic  refrigerator  or  magnet  1  e  heat  pump  ;  s.  tie  nagu.-toca  1  in  c  of  i  rot 

of  a  solid  material  to  e>t  ract  hear  f  rum  a  low- 1  eir.porat  a  re  source  ,>-a  :  f  ••  r  i  »  t  '■ 

a  higher  temperature*  « f  uk.  The  entropy  change  govern!  ■  g  t.  is  pro.  «••••*  i  -.n.  l  i  or. 

of  temperature  (T  )  and  magnetization  (M).  The  total  der  i  v  it  i  ve  of  t  be  enl  r  n  .•  i;,  a 

magnetic  material  is  given  by  Kq.  1. 


«■.  «■ .  i'.  i'.  i'.  e.v.  i'.  i<.  ii.  i'.  i‘n'.  i 


as 

dT  +  -  dM 

3m  T 


Substituting  some  generalized  thermodynamic  relationships  [Wark,  7]  into 


Eq.  1  results  in  the  following: 


Cb  3M 

dS  =  _  dT  +  -  dB 

T  8T  R 


where  S  is  entrophy,  B  is  magnetic  field,  Cb  is  heat  capacity  at  constant  magnetic 
field,  M  is  magnetization,  and  T  is  the  absolute  temperature.  Equation  2  can  be 
solved  if  the  magnetic  field-dependent  heat  capacity  and  the  equation  of  state  for 
the  ferromagnetic  material  are  known.  The  equation  of  state  for  a  ferromagnetic 
material  can  be  calculated  by  using  the  molecular-field  approximation  [2],  The 
field  dependent  heat  capacity,  on  the  other  hand,  is  a  quantity  that  must  be  mea¬ 
sured  because  of  the  complex  crystalline  structures  of  the  ferromagnetic  material. 

For  an  isothermal  entropy  change,  Eq.  2  reduces  to  the  following  form: 


dS  =  (  3M/3T )B  dB 


and  for  a  reversible  adiabatic  temperature  change,  Eq.  2  reduces  to  Eq.  A. 


dT  =  - 


T  3M 


Ck  3T 


In  addition,  a  paramagnetic  solid  that  follows  Curie's  law  will  result  in 
(3M/3T)g  being  negative,  and  thus  one  can  make  some  qualitative  observations. 
Equation  3  indicates  that  as  the  field  is  changed,  an  inverse  isothermal  entropy 
change  occurs.  This  experimental  characteristic  is  illustrated  in  Fig.  1  for  a 


Wscv 


typical  ferromagnetic  material.  Also,  Eq.  4  indicates  that  as  the  field  Is 
changed,  a  direct  adiabatic  temperature  change  occurs,  that  is,  an  increase  in  field 
increases  the  temperature,  while  a  decrease  in  field  decreases  temperature.  This 
phenomenon,  which  is  illustrated  in  Fig.  2  for  a  typical  ferromagnetic  material, 
is  known  as  the  magnetocaloric  effect.  It  is  the  principle  upon  which  the  devel¬ 
opment  of  a  magnetic  refrigerator  is  based. 
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TEMPERATURE 

Fig.  1.  Entropy-temperature  characteristic  for  a 
magnetocaloric  effect.  (B  is  magnetic 
field,  T  is  temperature). 

If,  then,  the  magnetocaloric  effect  is  applied  in  a  cyclic  process,  one  can 
obtain  magnetic  cycles  equivalent  to  the  familiar  thermodynamic  gas  cycles.  The 
most  familiar  cycle  and  the  basis  for  comparing  all  other  cycles  is  the  Carnot 
cycle.  This  cycle  consists  of  two  isentropic  and  two  isothermal  steps  which  result 
in  heat  being  pumped  from  a  low-temperature  source  to  a  high-temperature  sink. 
Similarly,  in  a  magnetic  Carnot  cycle,  the  heat  is  transferred  from  the  low-tempera¬ 
ture  source  into  the  high-temperature  sink  in  isothermal  steps.  The  Carnot  cycle 


is  illustrated  schematically  in  Fig.  3a.  For  the  magnetic  case  this  cycle  can  span 
only  about  half  the  adiabatic  temperature  change  that  the  ferromagnetic  material 
undergoes,  which  is  typically  only  about  8°C  at  the  material  s  Curie  point. 


TEMPERATURE 

Fig.  2.  Adiabatic  temperature  change  vs.  temperature  for  a  ferro¬ 
magnetic  material  having  a  raagnetocaloric  effect. 

Larger  temperature  spans  require  a  cycle  with  some  type  of  heat  exchange.  One 
such  cycle  is  the  Brayton  cycle.  Figure  3b  illustrates  this  cycle  and  compares  the 
use  to  the  magnetic  cycle.  The  gas  cycle  process  uses  two  isentropic  steps  and  two 
isobaric  steps,  while  the  magnetic  process  has  two  isentropic  steps  and  two  iso¬ 
field  (constant-field)  steps.  The  heat  exchange  process  takes  place  during  the  iso- 
baric  steps  in  the  gas  cycle  and  the  isofield  steps  in  the  magnetic  cycle.  The 
isentropic  magnetization  and  demagnetization  processes  in  the  magnetic  cycle  are 
equivalent  to  the  ideal  compression  and  expansion  processes,  respectively,  in  the 
gas  cycle.  In  a  gas  cycle  the  expansion  and  especially  the  compression  stages  of 
the  cycle  are  not  typically  isentropic,  but  the  magnetic  process  is  isentropic. 
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Fig.  3a.  Carnot  Cycle 
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Fig.  3b.  Brayton  Cycle. 


Thermodynamic  refrigeration;  Carnot  and  Brayton  cycles. 
(T  =  temperature,  S  =  entropy,  B  =  magnetic  field, 
and  P  =  pressure). 


The  isentropic  nature  of  the  process  is  the  basis  for  the  much  higher  efficiency 
possible  in  magnetic  refrigeration.  In  fact,  the  power  required  for  a  magnetic 
cycle  has  been  estimated  to  be  50-70%  of  Carnot,  as  compared  with  5  to  10%  of 
Carnot  for  a  gas  cycle.  Thus  a  great  efficiency  advantage  exists  for  a  magnetic 
refrigerator  if  efficient  heat  transfer  is  obtainable. 

In  any  cycle  comparison,  the  magnitude  of  the  refrigeration  cooling  that  takes 
place  in  each  cycle  is  also  important.  The  isentropic  temperature  change  for  a 
typical  magnetic  refrigerator  is  only  16 °C,  while  a  gas  cycle  refrigerator  expanding 
gas  from  15  to  1  atmosphere  near  room  temperature  can  result  in  an  isentropic  tem¬ 
perature  change  of  about  150°C.  The  refrigeration  capacity  of  these  two  cycles  is 
best  compared  by  using  the  isothermal  entropy  change.  In  a  comparison  based  on  a 
unit  of  mass,  the  gas  cycle  is  better  than  the  magnetic  cycle  by  a  factor  of  20  near 
room  temperature,  but  the  advantage  decreases  with  decreasing  temperature. 

If  the  two  cycles  are  compared  based  on  the  volume  they  require,  the  magnetic 
refrigerator  at  room  temperature  is  better  by  a  factor  of  50,  because  the  density  of 
a  ferromagnetic  solid  is  1000  times  greater  than  the  density  of  a  gas  at  standard 
conditions.  The  conclusion  would  be  that  ideally  the  magnetic  refrigerator  is  a 
more  compact  process.  However,  this  compactness  presents  the  problem  that  the 
heat  exchanger  must  operate  in  a  much  smaller  volume.  In  a  small  volume,  the 
required  high  heat  exchanger  effectiveness  is  much  more  difficult  to  obtain. 

MODELING  THE  MAGNETOCALORIC  EFFECT 

The  magnetic  refrigeration  process  utilizing  the  raagnetocaloric  effect  on  a 
ferromagnetic  material  was  investigated.  The  primary  focus  of  the  study  was  to 
obtain  an  understanding  of  heat  transfer  and  more  specifically  the  regeneration  pro¬ 
cess  in  a  magnetic  Brayton  cycle. 


THE  CONCEPT 


The  experimental  apparatus  used  to  study  a  regenerator  in  a  magnetic  refrig¬ 
erator  is  shown  schematically  in  Fig.  4.  The  magnetic  refrigerator  consisted  of 
five  major  parts:  (1)  the  high  field  system,  consisting  of  a  superconducting  Magnet 
surrounded  by  a  pool  of  liquid  helium  and  contained  in  a  vacuum-insulated  Dewar 
flask;  (2)  the  ferromagnetic  material,  gadolinium,  a  rare-earth  metal  having  a 
Curie  point  of  293  K;  (3)  the  displacer,  which  moves  the  fluid  through  the  ferro¬ 
magnetic  material  to  perform  the  heat  exchange;  (4)  the  warm  heat  exchanger,  which 
removes  the  heat  produced  in  the  adiabatic  magnetization  process;  and  (5)  the 
copper  regenerator,  which  isolates  the  cold  section  of  the  refrigerator  from  the 
warm  mixing  chamber  of  the  displacer. 


Fig.  4.  Magnetic  refrigeration  concept. 


The  configuration  used  for  the  active  magnetic  material  in  the  model  was 
strips  of  embossed  gadolinium  ribbon  wound  into  nearly  flat,  circular  sections 
(called  "pancakes")  and  then  stacked  together  so  as  to  fill  a  fiberglass  sleeve. 

This  configuration  allowed  maximum  contact  area  between  the  magnetic  material  and 
the  working  fluid,  and  could  in  fact  be  fabricated  when  the  model  was  verified 
experimentally  (see  also  Green  et  al.)  [8]. 

In  the  model,  the  magnetic  cycle  takes  place  as  follows:  (1)  The  magnet  is 
ramped  up  to  full  field,  causing  the  ferromagnetic  material  to  heat  up.  (This  step 
is  equivalent  to  the  compression  process  in  a  gas  cycle.)  (2)  The  displacer  moves 
down,  pushing  the  working  fluid,  nitrogen,  through  the  ferromagnetic  material  and 
removing  the  heat  of  the  magnetization  process.  This  heat  is  carried  to  the  warm 
heat  exchanger  and  carried  off  by  the  counterflowing  water  circuit.  (3)  The 
magnetic  field  is  removed,  cooling  the  gadolinium  through  an  adiabatic  demagnetiza¬ 
tion  process.  (4)  The  displacer  is  moved  up,  and  the  fluid  is  cooled  as  it 
passes  through  the  ferromagnetic  material  cooling  the  volume  between  the  cold  end 
of  the  gadolinium  and  the  copper  regenerator. 

THE  MATHEMATICAL  MODEL 

The  ferromagnetic  material,  gadolinium,  can  produce  a  significant  adiabatic 
temperature  change  when  subjected  to  a  high  magnetic  field.  An  adiabatic  tempera¬ 
ture  change  of  14°C  was  measured  by  Benford  and  Brown  [9]  near  its  Curie  temperature 
under  a  7-T  field  (Fig.  5).  The  magnetic  refrigerator  is  based  on  this  adiabatic 
temperature  change  and  the  thermal  exchange  mechanism  of  the  regeneration  process, 
which  also  occurs  in  the  ferromagnetic,  or  "active,"  material. 
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Fig.  5.  Magnetocaloric  characteristics  of  gadolinium. 

To  help  us  better  understand  this  process,  we  formulated  a  numerical  model  of 
the  magnetic  field  and  thermal  interactions.  The  magnetic  refrigerator  involves  four 
steps  per  cycle  to  obtain  the  cooling  effect  in  a  cyclic  process.  They  are: 

•  Magnetize  the  material 

•  Remove  the  heat  of  magnetization 

•  Demagnetize  the  material 

•  Cool  the  load 

The  magnetic  refrigeration  cycle  is  a  continuous  thermal  transient  process.  One 
could  visualize  this  transient  process  as  a  thermal  wave  that  passes  through  the 
active  material  every  half  cycle.  The  heat  of  magnetization  is  removed  and  the  load 
cooled  by  a  pressurized  gas  passed  through  the  magmetic  material  periodically  every 
half  cycle. 

A  2-dimensional  triangular  irregular  mesh  computer  program  (called  TRIM)  [10] 
was  used  to  calculate  the  magnetic  field  produced  by  the  superconducting  magnet 

indicated  in  Fig.  4. 
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Initial  calculations  assumed  no  ferromagnetic  material  in  the  bore  of  the  magnet. 
This  configuration  was  impractical,  however,  because  the  ferromagnetic  material 
(gadolinium)  had  certain  magnetic  properties.  The  permeability  of  gadolinium  was 
calculated  based  on  the  classical  molecular  field  theory  described  in  the  appendix. 
The  permeability  of  the  gadolinium  as  a  function  of  field  was  then  used  in  the  TRIM 
program  to  calculate  the  magnetic  field  in  the  gadolinium.  To  simulate  the  geometry 
of  the  gadolinium  a  series  of  0.0635-cm-thick.  concentric  rings  spaced  0.0635  cm 
apart  were  modeled.  The  field  distribution  with  and  without  and  concentric 
gadolinium  rings  in  the  bore  of  the  magnet  is  illustrated  in  Fig.  6. 


effect  on  the  thermal  performance  of  a  magnetic  refrigerator,  and  it  was  minimized 
in  the  experimental  apparatus  by  use  of  a  magnet  that  was  long  in  relation  to  the 
gadolinium  element. 


The  magnetic  refrigerator  performance,  in  which  the  gadolinium  undergoes  the 
adiabatic  temperature  change  and  performs  the  regeneration,  was  modeled  using  a 
finite  element  program  NASTRAN  111).  in  our  study  only  the  gadolinium,  the  cyclic 
fluid  motion,  and  the  stainless  steel  container  were  modeled.  The  model,  was  con¬ 
sidered  to  have  a  uniform  radial  temperature  gradient  and  thus  a  1-dimensional 
thermal  geometry.  We  also  assumed  that  the  temperature  entering  the  gadolinium 
from  the  warm  end  was  constant  (room  temperature)  and  that  the  temperature  entering 
the  gadolinium  from  the  cold  end  was  the  weighted  average  of  the  initial  and  final 
temperature  of  the  gas  exiting  from  the  gadolinium  during  previous  half  cycle. 

The  other  boundary  conditions  were  assumed  to  be  adiabatic.  The  partial  differential 
equation  for  the  gadolinium  and  the  stainless  steel  cylinder  is  as  follows: 


pm  =  Cm 
m  m 
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and  the  partial  differential  equation  for  the  fluid  is  given  by  Eq .  6: 
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In  addition,  the  following  assumptions  were  made: 


•  The  magnetic  field  in  the  gadolinium  calculated  with  TRIM  is  correct. 

•  The  heat  capacity  of  gadolinium  as  a  function  of  temperature  calculated  by 
Griffel  et  al.  (12]  is  correct. 


•  Thermal  conductivity  In  gadolinium  Is  not  a  function  of  temperature. 

•  The  adiabatic  demagnetization  temperature  change  as  measured  by  Benford  and 
Brown  [9]  is  accurate. 

•  The  working  fluid  obeys  the  Ideal  gas  law. 

•  The  heat  transfer  coefficient  determined  using  the  gap  regenerator  data  from 
Kays  and  London  [13]  is  accurate. 

•  The  contact  area  between  layers  of  gadolinium  is  about  one-tenth  of  the  cross 
sectional  area. 

•  The  ineffectiveness  and  heat  capacity  of  the  copper  regenerator  is  not  signi¬ 
ficant  for  the  isothermal  boundary  condition. 

The  finite  element  math  model  consisted  of  eight  elements  for  the  gadolinium, 
eight  elements  for  the  working  fluid  (nitrogen  gas)  and  eight  elements  for  the 
stainless  steel  container.  A  transient  thermal  analysis  was  performed  between  each 
magnetization  and  demagnetization  step  of  a  cycle,  assuming  a  uniform  field  of  5-T 
over  the  complete  length  of  the  gadolinium  material  during  the  magnetization  step. 

The  baseline  performance  was  calculated  based  on  the  assumptions  listed,  with  the 
results  plotted  in  Fig.  -7. 

For  the  second  scenario  modeled,  we  assumed  that  the  gadolinium  ribbon  was 
sliced  into  pancakes  thinner  than  the  original  3.8  cm,  and  an  equivalent  thermal 
resistance  was  calculated  based  on  0.32-cm  thick  pancakes  having  a  total  stack 
length  of  15.24  cm.  This  configuration  gave  0.07  the  thermal  conductance  compared 
with  the  original  baseline  case.  Although  the  improvement  was  small  for  this 
case,  a  larger  temperature  span  would  amplify  the  advantage  of  the  reduced  thermal 
conductance. 

In  the  third  scenario,  we  again  assumed  resistance  based  on  the  thinner 
gadolinium  ribbon  pancakes  and,  in  addition,  assumed  a  uniform  7-T  magnetic  field 
over  the  complete  length  of  the  gadolinium  material.  This  drastically  changed  the 
refrigeration  capacity  and  the  available  temperature  span  obtainable  using  gadolinium. 
In  fact,  a  60-K  temperature  span  was  calculated.  This  result  is  compared  to  the 


previous  two  cases  in  Fig.  7.  "he  third  scenario  for  cold  end  temperature  predic¬ 
tions  calculated  with  this  finite  element  model  agreed  closely  with  the  measurements 
made  by  Brown  [5]  in  his  experiments  using  a  7-T  field.  These  calculations  indicate 
the  significance  of  having  a  high  magnetic  field  throughout  the  length  of  the 
gadolinium. 


CYCLES 

Fig.  7.  Calculated  magnetic  refrigeration  cooling  characteristics. 

For  the  second  scenario  modeled,  the  Gd  ribbon  was  sliced 
into  thinner  sections,  producing  a  thermal  conductivity 
0.01  times  that  of  the  baseline  model;  for  the  third 
scenario,  a  constant  magnetic  field  of  7T  was  used  in 
addition  to  the  thinner  Gd  ribbon. 

To  span  even  larger  temperatures  in  a  single  stage  of  a  magnetic  refrigerator, 
several  active  magnetic  elements  could  be  stacked  on  top  of  each  other  in  order  of 
decreasing  Curie  temperature.  This  would  result  in  a  greater  average  adiabatic  tem¬ 
perature  change  in  the  demagnetization  process  over  a  much  wider  temperature  span. 
Three  materials  that  might  have  these  characteristics  are  gadolinium  (Gd)  [9], 
terbium  (Tb),  and  dysprosium  (Dy)  (15].  The  raagnetocaloric  effect  for  these  is  illus¬ 
trated  in  Fig.  8.  The  characteristics  shown  for  terbium  were  assumed  to  be  similar 
to  gadolinium  and  dysprosium  due  to  terbium's  similar  heat  capacity  characteristics. 
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Fig.  8.  Magnetocaloric  characteristics  of  three  rare-earth  metals. 

These  three  materials  were  assumed  to  have  the  properties  shown  in  Fig.  8. 

They  were  entered  into  the  finite  element  model  as  a  combination  of  several  stacked 
active  magnetic  elements  and  their  performance  calculated  assuming  a  constant  7T 
magnetic  field  over  the  full  length  of  the  active  magnetic  elements.  The  most 
important  characteristic  of  this  particular  configuration  was  the  heat  capacity  of 
the  rare-earth  metals  in  the  cool-down  process.  The  cool-down  time  was  greatly 
increased  by  the  additional  heat  capacity  of  the  rare-earth  metals,  which  had  a 
small  magnetocaloric  effect  at  the  higher  temperatures.  In  order  to  optimize  the 
small  magnetocaloric  effect  and  minimize  the  heat  capacities  of  the  inactive  rare- 
earth  metals  while  they  are  being  cooled  to  their  Curie  points,  half  the  volume  of 


EXPERIMENTAL  VERIFICATION  OF  THE  MODEL 


The  experimental  apparatus  needed  to  measure  the  thermal  performance  of  the 
magnetocaloric  effect  in  an  active  regenerator  was  designed  and  fabricated  at 
DTNSRDC.  This  prototype  refrigerator  was  designed  to  simulate  a  reciprocating  mag¬ 
netic  refrigerator  having  a  single  stage,  and  its  operation  was  similar  to  the  con¬ 
cept  shown  in  Fig.  4  and  described  above  under  "Modeling  the  Magnetocaloric  Effect. 
The  measurements  obtained  from  experimental  tests  on  the  refrigerator  were  compared 
to  results  of  the  mathematical  model  discussed  in  the  previous  section. 

The  experimental  apparatus  consisted  of  five  basic  components  as  follows: 

1.  Superconducting  magnet 

2.  Heat  exchanger 

3.  Passive  regenerator 

4.  Active  ferromagnetic  regenerator 

5.  Displacer 

The  superconducting  magnet,  fabricated  using  a  niobium-titanium  supercon¬ 
ductor  in  a  copper  matrix,  consisted  of  some  7122  turns.  This  magnet  was  designed 
in  a  solenoid  configuration  having  a  length  1-2/3  times  that  of  the  ferromagnetic 
material,  assuring  that  the  magnetic  field  in  the  ferromagnetic  material  would  be 
reasonably  uniform  and  have  a  magnitude  of  about  7  T.  The  magnet  was  25.4  cm  long. 
The  bore  diameter  was  8.31  cm  and  the  outside  diameter  12.13  cm.  The  solenoid  mag¬ 
net  had  an  inductance  of  1.51  H  and  it  could  be  ramped  up  to  225  A  with  a  12-V 
power  supply  in  approximately  28  <’econds. 

Figure  10  shows  the  superconducting  solenoid  magnet  attached  to  its  support 
structure  and  the  magnet  connected  to  a  pair  of  vapor-cooled  leads.  A  cryogenic 
Dewar  flask  with  a  warm  bore  of  5.46  cm  was  constructed  to  hold  the  magnet  and  pro¬ 
vide  the  vacuum  insulation  for  the  liquid  helium.  The  magnet's  Dewar  flask  had  two 


vapor-cooled  shields  at  20  and  80  K  to  keep  the  heat  load  to  the  helium  at  a  mini¬ 
mum.  A  steady-state  boil-off  rate  of  1  liter  of  liquid  helium  per  hour  was  mea¬ 
sured  as  the  heat  load  conditions  for  the  Dewar  flask.  Figure  11  shows  the 
calculated  and  measured  magnetic  field  distribution  along  the  axis  of  the  magnetic 
solenoid  with  the  magnet  current  at  225  A.  Superimposed  on  Fig.  11  is  a  represen¬ 
tation  of  the  active  ferromagnetic  regenerator  (gadolinium).  The  figure  shows  that 
the  field  distribution  across  the  gadolinium  was  not  completely  uniform.  The  field 
varied  about  12%  from  the  center  to  either  end  of  the  gadolinium. 


Fig.  10.  Superconducting  magnet  used  in  the  experimental 
magnetic  refrigerator. 
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Fig.  11.  Magnetic  field  distribution  along  the  axis 
of  the  superconducting  magnet. 

(TC#  ■  Thermocouple  number). 

Most  of  the  magnetic  refrigerator  apparatus  is  shown  schematically  in  Fig.  12. 
The  magnetic  material  was  centered  in  the  magnetic  field  by  placing  this  apparatus 
through  the  warm  bore  of  the  magnet's  Dewar  flask.  The  passive  regenerator  and  the 
active  ferromagnetic  regenerator  both  were  constructed  as  narrow  embossed  metal 
ribbon  wound  into  discrete  sections  ("pancakes")  of  material  and  stacked  into  a  G-10 
fiberglass  sleeve  (see  Fig.  13).  The  passive  regenerator  was  fabricated  with  copper, 
the  magnetic  regenerator  with  gadolinium.  The  gadolinium  ribbon  was  embossed  using 
a  rolling  mill;  the  process  is  described  in  detail  by  Lare  [16].  The  formed 
gadolinium  ribbon  was  sheared  into  0.317-cm  widths  and  lengths  of  it  were  wound 
into  the  pancake  sections.  The  winding  mechanism  is  described  by  Patton.* 

These  pancakes  could  then  be  placed  into  the  fiberglass  cylinder  as  illustrated 
in  Fig.  13. 

♦Patton,  W.G.,  "Winding  Mechanism  for  Winding  Ribbon-Gap  Regenerators,"  Patent 
Disclosure  Application  (1985). 
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Fig.  12.  Magnetic  refrigerator  apparatus. 


Fig.  13.  Gadolinium  ribbon  regenerator  section  in  the 
G-10  fiberglass  sleeve. 


The  passive  regenerator  used  a  copper  wire  that  also  was  formed,  wound  into  a 


pancake,  and  inserted  into  a  fiberglass  cylinder.  The  embossed  bumps  on  the  metal 
ribbons  provided  a  set  of  uniformly  spaced  flow  channels  separated  by  a  high  heat 
capacity  material  to  provide  the  regeneration.  The  copper  regenerator  had 
0.0178-cm  embossed  bumps  on  a  ribbon  0.0254-cm  thick.  This  provided  a  packing  factor 
of  65%  and  a  total  copper  mass  of  2.47  kg.  The  gadolinium,  in  contrast,  had 
0.0127-cm  embossed  bumps  on  a  ribbon  0.0127-cra  thick,  which  produced  a  packing  factor 
of  50%  and  a  total  mass  of  0.900  kg. 

A  finned  tube  heat  exchanger  (Fig.  14)  was  placed  at  the  warm  end  of  the 
refrigerator  to  remove  the  heat  of  the  magnetization  process  and  to  maintain  a  con¬ 
stant  tecqjerature  at  that  point.  Water  flowing  through  the  tube  intercepts  the  heat 
from  the  working  fluid  that  passes  through  the  fins  of  the  shell  portion  of  the  heat 
exchanger. 


Fig.  14.  Finned  tube  heat  exchanger. 

In  addition,  a  small  heater  capable  of  producing  about  20  W  of  power  was  placed 
in  the  cold  section  of  the  apparatus.  The  heater  was  shown  on  Fig.  12.  It  was 
intended  for  use  in  measuring  the  refrigerator's  performance. 


The  outlet  portion  of  the  finned  tube  heat  exchanger  and  the  copper  regenerator 
were  connected  to  the  displacer  cylinder  on  either  side  of  the  piston.  This  dis¬ 
placer  was  a  hydraulic  cylinder  with  a  12 . 7-cm-diameter  piston  and  a  22.86-cm 
stroke.  Coupled  directly  to  the  displacer  was  a  hydraulic  cylinder  with  a  5.08-cm 
bore  and  22.86-cm  stroke;  this  cylinder  drove  the  displacer.  A  Rexroth  hydraulic 
pump  with  a  four-way  proportional  valve  moved  the  displacer  back  and  forth  under 
electronic  control,  thus  moving  the  gaseous  working  fluid  (pressurized  N2  gas) 
back  and  forth  through  the  magnetic  refrigerator. 

Temperatures  were  measured  in  the  refrigerator  components  with  type  T  (copper- 
constantin)  thermocouples.  These  thermocouples  were  inserted  through  small  holes 
drilled  in  the  fiberglass  tube  and  into  the  component  material.  Number  36  thermo¬ 
couple  wire  was  used  to  minimize  heat  conduction  and  response  time.  The  thermocouple 
wire  was  placed  in  grooves  machined  in  the  outside  of  the  fiberglass  tube.  Figure  15 
shows  some  of  the  thermocouples  in  the  magnetic  material  leading  to  the  right  end  of 
the  assembly,  where  they  pass  through  the  pressure-tight  feed-through.  The  location 
of  the  thermocouples  in  the  magnetic  material  is  shown  in  Fig.  11,  superimposed  in 
relation  to  the  magnetic  field.  A  total  of  13  thermocouples  were  used  as  follows: 

•  3  in  the  passive  regenerator 

•  1  on  the  heater 

•  1  in  the  gas  space  between  the  regenerator  and  the  magnetic  material 

•  7  in  the  active  magnetic  regenerator 

•  1  in  the  gas  space  between  the  magnetic  material  and  the  heat  exchanger 

The  refrigerator  was  controlled  and  the  data  recorded  using  a  small  digital 

laboratory  computer  equipped  with  the  necessary  1/0  modules. 
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Fig.  15.  Magnetic  refrigerator  assembly 


The  magnet  current  was  controlled  by  the  output  of  a  digital-out  line  which 
operated  the  ramp-enable  circuit  in  the  magnet  exciter.  The  current  level  for  the 
magnet  was  fed  back  to  an  analog-in  port  of  the  computer  using  an  analogous  voltage 
from  the  ammeter  on  the  exciter. 

The  displacer  was  controlled  in  a  similar  manner.  The  computer  program  oper¬ 
ated  through  an  analog-out  line  to  an  electronically  controlled  hydraulic  valve  to 
cause  the  displacer  to  go  up,  stop,  or  go  down.  Displacer  position  was  sensed  with 
a  linear  voltage  differential  transformer  (LVDT)  attached  to  the  displacer  shaft 
and  fed  back  to  the  computer  through  an  analog-in  port.  Displacer  speed  was  also 
programmable. 

Temperatures,  magnet  current,  displacer  position  and  cycle  number  were  recorded 
on  floppy  disks  and  simultaneously  displayed  on  the  computer  view  terminal.  The 
control  and  data  acquisition  program  was  written  in  four  segments,  which  corresponded 
to  the  four  phases  of  the  magnetic  cycle.  Subroutines  were  also  invoked  to  display 


The  mass  flow  of  the  working  fluid,  N2  gas,  was  varied  as  a  function  of  temper¬ 
ature  to  maintain  equivalent  changes  in  the  heat  capacity  of  the  gas  and  the  magneti 
material.  The  mass  flow  was  calculated  based  on  a  linear  approximation  of  the 
adiabatic  magnetization  process  and  is  plotted  in  Fig.  16.  The  mass  was  varied  by 
assuming  the  gas  pressure  would  remain  the  same  and  the  displacer  stroke  would  be 
varied  as  a  function  of  gadolinium  cold  end  temperature.  This  displacer  character¬ 
istic  could  then  be  controlled  by  the  computer  program. 


Fig.  16.  Calculated  mass  flow  characteristics  based  on  a  linear 
approximation  of  the  adiabatic  magnetization  process. 


The  magnet  exciter  was  capable  of  ramping  the  magnet  up  or  down  from  0  to  220  A 
in  approximately  30  s  with  an  ll-V  potential.  This  magnetization  process,  plus  the 
displacer  movement  and  the  printing  of  the  data,  resulted  in  a  90-s  elapsed  time 
span  for  each  cycle  of  operation.  Thus,  a  200-cycle  cool-down  test  from  room  tem¬ 


perature  would  require  approximately  5  h  of  operation. 


After  several  preliminary  experiments  to  determine  the  adiabatic  temperature 
change  of  the  material  as  a  function  of  temperature  with  a  7-T  change  in  magnetic 
field,  the  experimental  magnetic  refrigerator  was  used  to  perform  two  cool-down 
experiments  starting  at  room  temperature  (300  K).  One  run  lasted  8  h,  the  other 
5  h.  During  each  test,  the  refrigerator  was  allowed  to  operate  until  the  supply  of 
liquid  helium  to  cool  the  superconducting  magnet  was  depleted.  Between  tests  the 
refrigerator  was  allowed  to  warm  up  throughout,  and  the  displacer  was  cycled  up 
and  down  with  no  magnetic  field  present  to  make  sure  that  all  parts  of  the  regener¬ 
ator  and  active  material  had  reached  thermal  equilibrium. 

After  the  second  test,  several  problems  occurred  that  prevented  any  further 
tests.  A  leak  developed  in  the  heat  exchanger,  which  caused  all  but  two  of  the 
thermocouples  in  the  magnetic  material  to  fail;  leaking  water  from  the  heat  exchanger 
caused  movement  and  created  a  large  pressure  drop  in  the  ferromagnetic  material. 
However,  a  significant  amount  of  data  was  obtained  and  is  reported. 

RESULTS 

In  the  preliminary  experiments,  the  thermocouples  implanted  in  the  gadolinium 
were  used  to  measure  the  adiabatic  magnetization  and  demagnetization  process. 

Figure  17  compares  the  measured  temperature  change  resulting  from  the  magnetization 
step  in  the  gadolinium  with  previously  reported  data  [9].  Our  measurements  showed 
a  25%  reduction  in  the  performance  of  the  magnetic  material. 

The  short  fall  in  adiabatic  temperature  change  in  this  experiment  as  compared  to 
Benford  and  Brown's  results  can  be  attributed  to  several  factors: 

1.  About  15%  of  the  difference  in  adiabatic  temperature  change  can  be  attributed 
to  the  heat  capacity  of  the  nitrogen  gas  surrounding  the  gadolinium  ribbon. 

2.  There  is  some  heat  transfer  in  the  axial  direction  of  the  gadolinium  mate- 


ADIABATIC  TEMPERATURE  CHANGE  (K) 


3.  The  G-10  fiberglass  Cube  surrounding  the  gadolinium  has  some  heat  capacity 


which  diminishes  the  adiabatic  temperature  rise. 

4.  The  thermocouple  preamplifiers  used  to  interface  with  the  computer  data 
acquisition  system  had  at  least  a  1 %  order  in  output. 

These  four  factors  also  act  to  reduce  the  reversibility  of  the  adiabatic  tem¬ 
perature  change. 
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Fig.  17.  Measured  adiabatic  magnetization  characteristics 
for  Gd  ribbon. 
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The  two  cool-down  experimental  runs  conducted  with  this  apparatus  were  started 
near  room  temperature  (300  K)  and  cooled  the  cold  portion  of  the  refrigerator  down 


to  a  low  temperature  of  246  K.  Cool-down  profile  of  these  two  runs  is  shown  in 
Fig.  18.  These  two  cool-down  characteristics  showed  very  similar  performance;  the 
slight  variation  resulted  from  the  difference  in  the  mass  flow  of  nitrogen  gas  for 
each  test.  Run  2  was  conducted  under  more  nearly  optimal  mass  flow  conditions. 

In  one  of  the  early  preliminary  runs,  before  the  thermocouples  broke,  the  tem¬ 
perature  distribution  through  the  gadolinium  was  measured.  Each  set  of  measurements 
indicated  the  temperature  profile  after  each  significant  event  in  the  cycle. 

Figure  19  shows  these  characteristics  of  the  temperature  distribution  in  magnetic 
material.  In  addition,  the  figure  also  shows  the  variation  of  the  magnetocaloric 
effect  in  the  material  as  a  function  of  temperature  and  field.  This  smaller  mag¬ 
netocaloric  effect  at  the  lower  temperatures  could  be  improved  by  using  other  rare- 
earth  metals  with  decreasing  Curie  temperatures.  This  possibility  was  discussed  in 
the  section  on  the  mathematical  model  (see  also  Fig.  8). 
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Fig.  18.  Measured  cool-down  profile  of 
the  magnetic  refrigerator. 


The  average  measured  adiabatic  temperature  change  in  the  magnetic  material  was 


about  7°C.  The  lowest  temperature  obtained  with  the  experimental  apparatus  was 


246  K,  a  50°C  temperature  difference  across  the  active  regenerator.  These  measure¬ 


ments  compare  well  with  the  results  described  earlier  on  the  mathematical  model 


and  results  described  by  Brown  [5].  Although  the  math  model  did  not  include  the 


passive  regenerator  in  the  refrigerator  and  its  associated  capacity,  the  final  tem¬ 


perature  obtained  in  the  experimental  apparatus  was  close  to  that  calculated  in  the 
numerical  model. 
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Fig.  19.  Measured  temperature  distribution  of  the  Gd 
active  magnetic  regenerator. 


CONCLUSIONS 


The  experimental  magnetic  refrigerator  produced  a  50°C  temperature  span  with 


only  a  7°C  adiabatic  temperature  change  in  the  magnetic  material.  The  regeneration 


process  in  the  magnetic  material  produced  a  factor  of  7  increase  in  the  temperature 


span  for  a  single  ferromagnetic  material.  Without  a  doubt,  the  active  regeneration 


concept  used  in  conjucntion  with  the  embossed  ribbon  configuration  ^an  successfully 
produce  magnetocaloric  refrigeration  at  temperatures  higher  than  4  K. 


Further,  if  several  rare-earth  metals,  having  significant  magnetocaloric  effects 
at  decreasing  temperatures,  were  combined  as  the  active  magnetic  element,  the  possi¬ 
bility  exists  for  producing  a  single-stage  magnetic  refrigerator  that  could  operate 
over  a  large  temperature  span.  Such  a  refrigerator,  if  developed,  offers  the  possi¬ 
bility  of  a  compact,  highly  efficient,  reliable,  shock  resistant,  closed  cycle 
refrigerator  for  cooling  or  liquefying  gases  from  room  temperature  down  to  liquid 
helium  temperature  (4  K)  in  only  a  few  stages. 
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APPENDIX 


CLASSICAL  MOLECULAR  FIELD  THEORY 

The  strong  interaction  in  a  ferromagnetic  material  that  aligns  the  atomic 
dipoles  in  a  paraLlel  manner  may  be  represented  by  some  internal  field  HM.  The 
thermal  agitation  of  the  atoms  in  the  ferromagnetic  lattice  opposes  the  orientation 
of  the  field.  The  Curie  temperature  is  defined  as  the  point  at  which  the  thermal 
agitation  is  of  a  sufficient  magnitude  to  counteract  the  spontaneous  magnetization. 

Weiss  [17]  developed  a  phenomenological  model  for  ferromagnetism  by  assuming 
that  H>q  was  proportional  to  the  magnetization  M  (i.e.,  spontaneous  magnetization). 

HM  =  KM  (A .  1 ) 

The  symbol  K  is  termed  the  molecular  field  constant  (or  mean  field  coupling  con¬ 
stant).  The  field  HM  is  called  the  molecular  field  constant.  However,  the  field 
is  sometimes  termed  the  Weiss  or  exchange  field.  Ln  the  presence  of  an  arbitrary 

applied  field  H  the  actual  field  acting  on  a  particular  dipole  H’p  is 

HT  =  H  +  KM  (A. 2) 

The  demagnetizing  and  Lorentz  (dipole-dipole)  fields  are  ignored  because  they  are 
small  when  compared  to  the  molecular  field. 

SPONTANEOUS  MAGNETIZATION  REGION 

We  assume  that  the  atoms  of  the  solid  have  a  total  angular  momentum  quantum 

number  J  and  have  N  atoms  per  unit  volume.  Then  the  magnetization  of  the  ferro¬ 

magnetic  material  is  given,  in  Gaussian  units,  by:  [18] 


M  =  NgUgJ  Bj (x) 


(A. 3) 
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(A. 5) 


The  symbols  in  the  above  equations  have  the  following  meaning.  (The  numerical 
values  for  the  atoms  are  for  gadolinium). 

k  =  Boltzmann's  constant  (1.381  x  10“^  erg/K) 

J  =  total  angular  momentum 

K  =  mean  field  coupling  constant  (calculated) 


PR  =  bohr  magneton 


0.927  x  10  erg/gauss 


g  =  Lande  factor  (dimensionless  quantity) 

T  =  temperature  in  Kelvin 
B  =  magnetic  induction  (gauss) 

H  =  magnetic  field  (oersteds) 

M  =  magnetization  (magnetic  moment/cm^) 

Nq  =  Avogadro's  number  =  6.0225  x  102^  mole  ^ 

N  =  number  density  of  dipoles  (0.2596  x  10^3  Cm-^ ) 


T^  =  ferromagnetic  Curie  point 
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The  mean  field  constant  K  may  be  calculated  as  a  function  of  Tq,  the  Curie  tem- 


perture.  Above  Tq  no  spontaneous  M  exists  so  that  Eqs.  A. 3  and  A. A  have  only  a 
trivial  solution  for  M. 


NgzuB  J  ( J+l )K  =  3kTc 


(A. 6) 


Considering  the  above  theory,  the  equation  for  the  magnetization  M  of  gadolinium 

may  be  expressed  in  detail  as 
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Dividing  through  by  M  the  equation  is 


(A. 7) 
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(A. 8) 


For  a  particular  field  H  the  magnetization  of  the  gadolinium  may  be  calculated 
on  the  compute r  by  a  variational  procedure  which  picks  an  M  so  that  Eq.  A. 8  equals 
zero.  Once  a  particular  M  is  determined  for  a  chosen  H,  B  may  be  determined  from 
the  following  equation: 


B  =  H  +  AtiM 


(A. 9) 
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